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FOREWORD 


A study to determine the fesslbility of eccommodeting the R&D 
requirements of electroepltaxlal crystal growth in the Space Shuttle has 
been performed. The specific accomodations facilities studied were the 
Orbiter Middeck stowage lockers, Materials Experiment Assembly (MEA), and 
Get Away Special (GAS) Cans. The effort has encompassed development of 
guidelines and assumptions necessary to quantify and characterize elements 
of the electroepltaxlal process, conceptual design of a Gallium Arsenide 
Crystal Growth Facility, and assessment of facility requirements versus 
Orbiter Middeck, MEA, and GAS Can capability. 

The results of this study Indicates that the MEA can bast accommodate 
the R&D furnace facility. The Middeck area, though marginally suitable, 
has energy, heat rejection, volume, and safety concerns. The GAS Can 
program ground rules prohibit its use f^r this application. 

This study was performed by Teledyne Brown Engineering under 
Contract NAS8- 34743. 
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ELECTROEPIT AXIAL GROWTH OF GaAs 


Introduction: 

Crystal growth processes are phase transformations which are driven by differences 
in chemical potential between the nutrient and crystallizing phases. An important step in 
crystal growth is the transport of material from the bulk of the fluid to the growing 
crystal surface. In the case of melt growth this is achieved by temperature gradient 
which then dissipates the latent heat of fusion generated by the solidification process. 
Such temperature gradients in turn generate thermal convection in the melt. Both natural 
and forced convective flows provide very important soutces of enhancing mass transfer 
required for the growth of a crystal. 

One of the techniques used to achieve forced convection is to apply an electric 
field. This field may change the chemical potential of a species by either polarizing it 
and/or exerting a force on it equal to the product of the ionic charge of the species and 
the local electric field strength. 

The technique of electromigration, i.e., electric field induced forced convection, 
can be used to grow semiconductor material and other compounds from solution by passing 
electric current through the growth interface while the temperature of the system is 
maintained constant. Current-controlled electromigration, referred to as electroepitaxy, 
has been successfully applied to grow epitaxial layers of various semiconductors and 
garnets. 

Process Physics : 

An electric field can change the chemical potential of a component i by either 
polarizing it and/or exerting a force on it proportional to the ionic charge of i thereby 
causing a differential migration, i.e., diffusion. If ionized species are involved, 
electromigration due to field-charge interaction will typically dominate over the 
polarization term. The mass transfer under such conditions for the cations and anions are 
respectively. 

m c = +0+ c + FE 
m a = -B_ c_ FE 

Here, 0 is equivalent conductance, F is the Faraday, and E is the electric field vector. 
Thus the cations move in the direction of the electric field vector. In general, ions move 
at different rates and tht condition of electroneutrality must be 
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satisfied. This has several important consequences. Application of an electric field to an 
initially homogeneous solution causes partial separation of the ions. On the other hand, 
diffusion in a concentration gradient generates an ‘electric field. Diffusion of one ion into 
a homogeneous mixture of others can also cause their preferential migration. 
Crystallization in electrolytes causti electric fields to develop because o' different ionic 
diffusion rates, different tendencies towards incorporation of the different ions by the 
growing crystals, and differences in the absorption of anions and cations on the crystal 
surface. 

Mixtures of liquid metals are also found to exhibit electromigration. Here, however, 
the mechanism of relative movement is not as clear as it is for aqueous solutions. One 
can imagine the metal "ions" to be influenced by the field even though they are immersed 
in a cloud of electrons. However, recent work suggests that electromigration in liquid 
alloys is predominantly governed by a momentum exchange arising from nonelastic 
scattering of the electrons by the ions (1-3). 

From the experimental data it can be inferred that the electric current, rather than 
electric field, is responsible for elect* omigration in metals. Hence, the mass transfer flux 
equation should contain electric current terms as the driving force rather than electric 
field. Expe. : mentally, it has proven difficult to separate the effects of electromigration 
from the increased stirring generated by the Joule heating and magnetohydrodynamic 
forces caused by the high current densities (4). 

The growth rate during current-controlled electroepitaxy of GaAs can be explained 
on the basis of Verhoeven's Theory (5) for electromigration in metallic melt solutions. 
According to Verhoeven, the use of the differential mobility (y^) is only rigorously 
correct in dilute solutions. For more general case, the velocity of the solute is the 
resultant of the velocities due to diffusion, electromigration, and bulk movement of the 
interface. Under this condition the flux of the sol te towards the interface is given by 

* * - c a R * PM 5 b D ~3z~ + c a c bVd 2E 

where c a and Cfc, are the concentrations of the two components of the alloy, R is the rate 
of solidification, D is the mutual diffusion coefficient, p ^ is the molar density, Vb is the 
partial molar volume of the solvent component of the alloy, and E is the electric-field 
intensity. As an approximation, if the coefficients of c a are assumed constant, then the 
effective distribution coefficient can be calculated. 
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If the movement of the growth interface and the contribution due to diffusion can be 
neglected for GaAs electroepitaxy, then the As flux towards the interface can be 
approximated by 


“ c As c Ga^Ga y d 2E 

From a thermodynamic point of view, the key to the GaAs electroepitaxy is that the 
passage of electric current disturbs the equilibrium concentration of the gallium-arsenic 
solution by inducing the deposition of GaAs on to the GaAs substrate. Since the solution 
is saturated with arsenic at the equilibrium operating temperature, the removal of GaAs 
in the stoichiometric ratio will leave the solution arsenic deficient. To counterbalance, 
arsenic must be added to the solution. This is achieved through the desolution of a GaAs 
source material. 

The thermodynamic picture for the system is based on the binary phase diagram as 
depicted in Figure 1. The major, rather unique feature of the phase diagram is that a 
single liquidus curve covers nearly the entire composition range. The solid and liquid 
phases of the phase diagram down to quite low te peratures clarify the need and 
usefulness of the epitaxial growth technique. Moreover, the vapor phase plays a 
significant role in determining the optimum growth technique. For GaAs at its melting 
point, the arsenic pressure is of the order of 10 atm. On the Ga-rich portion of the 
liquidus at around 900 °C, the As pressure is < 1(H* atm as shown in Figure 2. This means 
that electroepitaxy can be carried out in a l-atm apparatus. 

Stoichiometric GaAs is formed as long as the rate of collision between absorbed 
gallium and impinging arsenic is equal to the arrival rate of gallium. Figure 3 shows 
schematically the deposition rate as a function of the ratio of arsenic and gallium arrival 
rates, for a fixed gallium arrival rate. Below RAs/^Ga = g row th is rich in 

gallium. Stoichiometric growth occurs at and beyond this ratio with a constant deposition 
rate until, in principle, arsenic begins to condense. 
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FIGURE 2. EQUILIBRIUM PARTIAL PRESSURES OF THE Ga-As SYST 
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FIGURE 3. DEPOSITION RATE VS. ARRIVAL RATE OF As 2 
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Present Laboratory Technique: 

The apparatus used in the laboratory studies is shown schematically in Figure 4. The 
growth cell is composed of two graphite segments each having a stainless steel electrode 
threaded into it (6). The two segments are electrically isolated with boron nitride, and 
current can flow through the system only when the GaAs substrate is brought in contact 
with the Ga solution with As in it. Substrates of 1.2 x 1.2 cm area, having thickness 
between 200 to 300 urn, with (100) crystallographic orientation were employed. A uniform 
electrical contact between the substrate and graphite was established with a 150-um- 
thic'< Ga layer. The growth cell i3 placed in a cylindrical furnace with goid reflectors for 
operational temperature control. 

In a typical growth experiment the cell is brought to the desired temperature and 
the galiium solution is equilibrated using a dummy GaAs substrate for a period ranging 
from 2 to 5 hr, depending on the quantity of the solution. The real substrate is then 
brought in contact with the solution, and electrical current is passed through the solution 
keeping the substrate positive with respect to the solution. The isolation of the real 
substrate during heat up and equilibration period was necessary to protect the substrate 
surface from damage. 

Experimental Results Important to Space Processing ; 

The experimented results pertinent to the space experiment facility design and 
development are presented below. In this context, the dependence of growth rate on 
current density, process temperature, and solution depths and the adverse effects of Joule 
heating and growth current fluctuations are discussed. The data on the dependence of 
growth rate on current density, etc., are needed for facility sizing and design, whereas the 
information on the effects of Joule heating and growth current fluctuations are necr lary 
for the facility development so that proper measures can be taken to avoid/minindze 
these detrimental effects. 

Effect of Current Density on Growth Rate - The growth was found to be increasing 
linearly with current density up to 60 A/cm 2 (7). This is shown in Figure 5 for a 
temperature of 900 °C and solution depth of 0.7 cm. The growth rate for Si-doped 
solution was found to be less than undoped solution at higher current density. 

Epilayers with the most uniform thickness were obtained at current densities of 50- 
60 A/cm 2 . At these current densities, the random fluctuations associated with 
nonuniform current density were counter-balanced by convective flow in solution due to 
localized Joule heating. 
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FIGURE it. LABORATORY FACILITY 
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FIGURE 5. EFFECT OF CURRENT DENSITY ON GROWTH RATE 
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Effect >f Temperature on Growth Rate - Figure 6 shows the dependence of the 
growth rate on temperature under a constant current density of 10 A/cm* (7). Below 
about 825 °C it is found that desolution of substrate taxes place. At this temperature, 
growth occurs only upon reversing die current polarity. This phenomenon is also true for 
the Si-doped solution. The increase in growth rate from 1.6 pm/min at 930 °C to 2.3 
ym/min at 973 °C may be an indication of better growth rate at higher temperatures. 

The dopant (Si) distribution coefficient increased with increasing temperature and 
increasing current density 

Effect of Joule Heating - Temperature measurements, after thermal equilibrium of 
the growth ceil, were carried out along a horizontal line extending to both sides of the 
growth cell (6). The results in the absence and presence of 20 amps of current flow are 
shown in Figure 7. 

Joule heating in electroepitaxy establishes a new thermal equilibrium characterized 
by steeper thermal gradients. In the experimental setup of Reference 6 the thermal 
gradients, after thermal equilibrium, are found to be a function of the electric power. 
This is shown in Figure 8. It has been found that Joule heating can play a significant role 
and introduce defect in the grown single crystal material. Localized Joule heating due to 
high contact resistance between cell and current electrodes can magnify the problem. 

Effect of Solution Depth or. Growth Rate - The effect of variation of the solution 
depth on growth rate is very small at smaller depths. This effect levels out for larger 
solution depths. This is shown in Figure 9 for a constant current density of 30 A/cm* and 
a temperature of 900 °C. 

Effect of Growth Current Fluctuation - Nonuniform current density at the growth 
interface will introduce random fluctuations in the epilayer thickness. Apart from source 
current fluctuation, the growth current fluctuations can arise from random defects in the 
eiectric contact to the substrate or fluctuations in the resistance across the substrate 
caused by fluctuations in the depth of substrate dissolution by the Ga contact layer. 

Gravity Related Crystal Growth Dynamics : 

As mentioned earlier, convective flows provide very important sources of enhancing 
material transfer from the fluid to the growing interface of the crystal. These convective 
flows can either be natural or forced. There is a multitude of driving mechanisms for 
natural convection including buoyancy and surface tension. In buoyancy driven convection 
fluid motions are caused entirely by the action of gravity field on density gradients in the 
fluid arising from temperature and/or concentration gradients. 
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FIGURE 6. EFFECT OF TEMPERATURE ON GROWTH RATE 
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Temperature Gradient Across Solution (°C/cm) 



FIGURE 8. TEMPERATURE GRADIENT ACROSS 
SOLUTION VS. ELECTRIC POWER 
AT 902°C 
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FIGURE 9. EFFECT OF SOLUTION DEPTH ON GROWTH RATE 
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Thermal Convection - The orientation of the heating and/or cooling surfaces with 
respect to the gravitational direction is very important in thermal convection flows. The 
nature of the shear flows generated by both gravity-parallel and -non parallel heat flows is 
an important aspect of crystal growth configurations with regard to thermal convection. 
The thermally driven convective forces can be expressed in a dimensionless form of the 
Navier-Stokes flow equations. For heat flow parallel to the gravity field, g, the 
temperature difference AT (t , is expressed as a Raleigh number 



’ AT » 

<n 


where x is the thermal expansion coefficient, k is the thermal diffusivity, n is the 
kinematic viscosity, and is the convective length. 

For a gravity-normal temperature difference, ATj_, the driving force is expressed as 
the dimensionless Grashof number 


, 0Oi x *AT x 

N Gr * ii* 

As the driving force for thermal convection increases the convection cells develop 
instabilities that result in time-dependent temperature and velocity behavior. Theoretical 
results on stable and unstable flows for NR a /Nc r =l are given in -the regime diagram 
depicted in Figure 10 in terms of the Grashof number and the Rayleigh number. 

There are many factors which influence the critical Rayleigh or Grashof number at 
which instability occurs. In general terms, these factors include thermal properties, 
variations in melt properties with temperature, orientation with respect to gravity field, 
container geometry, and presence of other body forces. 

Chemical Convection - Density gradients due to compositional differences are 
important for solution growth environment. Solutal convection in isothermal systems is 
most easily studied by observing the dissolution of solids and measuring the associated 
mass transport rates. The diffusion boundary layer at a growth interface builds up by 
segregration processes and, at the point of gravitational instability, releases a blob of 
solute. The solutal release period can be computed from the solutal Rayleigh number and 
is predicted to vary with growth rate. 
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FIGURE 10. STABILITY DIAGRAM FOR CONVECTIVE FLOW 
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The coexistence of density gradients due to both temperature and composition 
(thermosolutal convection) is of more interest in crystai growth. 

Surface Tension Driven Convection - Surface-tension-driven convection is of 
particular interest in connection with materials processing efforts in the low-gravity 
environment of space. On Earth, surface-tension is important only for small-scale 
configurations, like thin-films, droplets, etc. Surface-tension depends on temperature, 
composition, and electric field. Horizontal temperature gradient, i.e., AT parallel to the 
free surface can sustain convective flow even in the absence of gravity field. The driving 
force for the flow is described by a dimensionless quantity called the Marangoni number 

N . . J2_ 

Ma <np dT 

where p is the fluid density, h is the layer height, and a is the surface tension. The 
critical NM a on Earth for a free surface with no heat flow is There is a tight 

coupling between Marangoni flow and buoyancy-driven flow such that the thermocapillary 
convection can be expressed as 

V M Ra + % /N 5a ■ 1 

where N& a and N^ a are the values at marginal stability for the two agencies that 
reinforce one another. Thus a liquid stable to buoyant forces can still be destabilized by 
surface tension forces. 


17 



References: 


1) S. G. Epstein and A. Paskin, Phys. Letts. 24A, 309 (1967) 

2) S. G. Epstein, Advances in Phys. U>» 323 (1967) 

3) P. P. Kuz'menko, Ye. I. Khar'kov, and V. I. Lozovoy, Phys. Metals and 

Metallog. 21, 94 (1966) 

4) 3. D. Verhoeven, Trans. Quart. Am. Soc. Metals 55, 866 (1967) 

5) 3. D. Verhoeven, Trans. Met. Soc. AIME 233 , 1156 (1965) 

6) L. 3astrzebski, Y. Imamura, and H. C. Gatos, 3. Electrochem. Soc. 125, 1140 

(1978) 

7) L. 3astrzebski and H. C. Gatos, 3. Cryst. Growth 42, 309 (1977). 


18 



